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Summary 

Calculations of inclusive inelastic scattering distributions for heavy-ion collisions are con- 
sidered within the high-energy optical model. Using ground-state sum rules, the inclusive 
projectile and complete projectile-target inelastic angular distributions are treated in both 
independent particle and correlated nuclear models. Comparisons between the models intro- 
duced are made for 4 He particles colliding with 4 He, 12 C, and ie O targets and for protons 
colliding with 16 0 targets. Results indicate that correlations contribute significantly, at small 
momentum transfers, to the inelastic sum. Correlation effects are hidden, however, when total 
scattering distributions are considered because of the dominance of clastic scattering at small 
momentum transfers. 

Introduction 

Research efforts that consider the effects of high-energy nuclei, such as cosmic rays, on 
physical systems require a large nuclear cross-section data base as input. The importance 
of the accuracy of this data base for space-radiation transport calculations has recently been 
discussed by Townsend and Wilson (refs. 1 and 2). The energy dependence and normalization 
of nuclear fragmentation parameters were shown to have a large effect on the prediction of 
particle differential flux and absorbed dose in tissue and aluminum shielding. An accurate 
fragmentation data base is therefore essential in order to assess the effects on astronauts from 
these radiations. 

For high-energy, charged-particle transport, the straight ahead approximation is seen to 
be accurate (refs. 3 qind 4), and the fragmentation inputs reduce to energy distributions for 
the scattered primaries and projectile/target fragments. The large number of nuclei and the 
extended energy range of importance in these transport studies require theoretical predictions 
of these distributions that are accurate and comprehensive. The high-energy optical model, as 
derived from Watson’s multiple-scattering series (refs. 5 and 6) or the Glauber approximation 
(refs. 7 and 8), gives reliable predictions for both total and absorption cross sections (ref. 9) 
and includes the energy dependence explicitly in the optical potential. In this report, we 
extend that work by evaluating angular distributions for inelastically scattered primaries using 
the high-energy optical model. We also consider the effects of nuclear correlations on these 
distributions. 

The Glauber model has been used to study inclusive scattering for “elementary” projectiles 
(refs. 8 and 10) and for heavy-ion scattering in the rigid-projectile approximation (refs. 11 
and 12). Semiclassical cascade equations have also been derived with the Glauber formalism 
(refs. 13-15). We generalize the Glauber result for inclusive inelastic scattering in the 
independent particle model to the heavy-ion case. The eikonal-coupled channel (ECC) 
amplitude (refs. 5 and 6) can be considered a matrix representation of the Glauber amplitude 
(refs. 16 and 17) when correct or equivalent kinematics are assumed. Using the ECC, the 
inclusion of correlated nuclear basis functions for evaluation of inelastic sums is handled in a 
straightforward manner. Correlation effects will contribute to leading order in the inelastic 
sum rule and therefore may become important for small momentum transfers. We consider a 
nonperturbative approach for including the effects of two-particle correlations in the inclusive 
sum. Comparisons between the models introduced will be made for protons (p) and alphas (a?) 
scattering on 4 He, 12 C, and ie O targets at several energies. 

Inclusive Inelastic Scattering 

We consider nucleus-nucleus scattering at high energies for the case where an inclusive 
measurement of the projectile final state is made, 


P + T-+P + X 


( 1 ) 



with P and T denoting the projectile and target, respectively, and X denoting some final state 
of the target that is not measured. In equation (1) the projectile scatters elastically, and meson 
production is not considered. In the overall center-of-mass (CM) frame, with the projectile and 
target states denoted by jn P > and | v T >, respectively, the angular distribution for equation (1) 
is found by summing (TOT) the nuclear-scattering operator over all final states of the target, 


do p \ 

/ TOT 


= | < i^tOp\]{^)\OtOp 

up 


> 


2 


(2) 


where / is the scattering operator and q is the momentum transfer to the projectile defined by 

q = k — k p W 

In equation (3) k and kp are the initial and final projectile wave vectors, respectively. In 
equation (2) the phase space is approximated by a two-body phase space that is expected to be 
accurate at high energies. Equation (2) can be separated into elastic (EL) and inelastic (IN) 
contributions given by 

—) = | < OpO T |/(q)|OpO T > | 2 ( 4 ) 

'el 


and 


dnj 

da p 
dft 


\ = 5Z I < Opn T |/(q)|OpOp 

/ IN 


> 


( 5 ) 


respectively. The summation in equation (5) includes all excited states, bound and continuum, 
of the target. This infinite summation can be reduced to a single matrix element through the 
use of closure on the target states: 


y \vp >< i>t\ — i — i 0t >< °t\ 

I/7' 


( 6 ) 


Inserting equation (6) into equation (5), gives 

da p 


d<7 P \ __ da p \ _ do \ 

^ / IN TOT ^ ' 


( 7 ) 


EL 


where 


do 


dQ 


=< 0 T | < 0 P |/(q)|0p X Op|/ + (q)|Op > 1 0 T > 


( 8 ) 


TOT 


The great advantage of equation (7) over equation (5) is that only the ground-state wave 
function of the target is needed. 

A second reaction that we consider is complete inelastic scattering 


P + T -* X + Y 


( 9 ) 


where the projectile and target are both left in excited states (denoted by PT). The angular 
distribution for equation (9) is given by 


da 


PT' 




| = ^ | < t'r^p|/(q)l°P C) T > I 2 

IN np^O 


( 10 ) 
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which is written, using closure on both the target and projectile states, as 


da PT \ _ da\ da\ da p \ da T \ 

d / IN dSl / tot + dfl ) £l dQ J T()T dU J TOT 

where 

^) T0T =< °pOt\ (l/(q)l 2 ) I o P o T > 

Equation (11) may be written as 



( 11 ) 


( 12 ) 


(13) 


The distributions considered above are evaluated when models for the nuclear-scattering 
operator and ground-state wave functions are introduced. 


Glauber Independent Particle Model 

The Glauber scattering operator (ref. 8) is defined in terms of the nucleus-nucleus profile 
function as 

A lie f ^ 

/(q) = — y exp(iq - b) T(b) d 2 b (14) 

where b is the impact parameter, and the profile function is given by 

Ap,Ap 

r (b) = 1 - P [1 -r aj -(b-s Q --Sj)] (15) 

Otj 

where T a j is the two-body profile function, a and j label the projectile and target constituents, 
respectively, and s is the projection of the internal nuclear coordinate onto the impact parameter 
plane. 

In the independent particle model (IPM), the nuclear transition density is approximated by 
the product of single-particle densities, 


Ap 

Pv t Ot( t 1 ■ 1 1 PvpOp ) ( 16 ) 

3 = 1 


with 


J' Pi^pOp (r)c?r S t /j’Op 


(17) 


The ground-state single-particle densities poo of the projectile and target are denoted by pp 
and pp, respectively. 

We now consider the evaluation of the distributions introduced in the preceding section 
using the Glauber IPM. For elastic scattering, we have from equation (14), 


^ 1 Jc f r \ ^ 

< O pOp\f (q)|0 pOj >= — J d 2 b < OpO T \r(b)\OpO T > exp(fq • b) (18) 


Then, we find from equation (4), 


da\ 
dfi) el 


^1 J exp [tq • (b - b')] d 2 b d 2 b' |l - [J t 1 - M *j( b )] } j 1 " II [ 2 ~ M aj( h> )} } 

^ otj } ^ Otj * 

(19) 
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where 


( 20 ) 


M a j( b) = J dr a dTj pp(r a ) Pr( r j) r aj( b “ s <* ~ s j) 


For inelastic scattering, we insert equations (14) and (15) into equation (8) to find 

= |**|2 [ fl b fly exp [iq.(b-b')] < Or I < Op\f(b)\Op 

d n J TOT 2?r d 

x < O P [f + (b')|0p > 1 0 T > 

After completing some algebraic steps, we find 

da F ^ 


> 


( 21 ) 




TOT 


: ft 2 [d 2 b dh'ex p [iq ■ (b - b')] (l ~ n [1 - *M b )] - II I 1 - M aj( h ')} 

2n J l aj aj 

+ P j) — M a j{h) — M£j( b ) + fi a j(b, b ) | 

aj ' 


( 22 ) 


where 


! aj (b,b') = J dTj pj{Tj) J dr a pp(r a )T a j(b-s a ~Sj)xY[ J dr a'Pp( r a') T^(b S a i 


Q 

Using equations (19) and (22) in equation (7) gives the inclusive inelastic distribution 
da p _\ I ik 12 

r, IN 


s j) 

(23) 


as 


dfi 


^-| j d 2 b d 2 6 7 exp [iq ■ (b - b 7 )] | JJ £l — Maj(b) - M+j(b ) + fi Q j(b,b )j 


n [i - i 1 - M U h 

aj 


■')]) 


(24) 


Thus, the form for the inclusive nucleus-nucleus distribution, as given by equation (24), is 
identical to the nucleon-nucleus case found in reference 8 with the exception of the form for 
Sl aj given by equation (23). Similar to the nucleon-nucleus case (ref. 8), equation (24) can be 
approximated in the large mass number and large momentum transfer limits by 


dQ 

where 


and 


A = ^-1 J d 2 b d 2 b f exp [iq • (b — b 7 )] exp ji |x(b) — X + (b )] } { ex P [^(b,b )] l} 
/IN for. 


ix( b ) = 


(25) 

(26) 


n(b,W) = Y, n aj(b,b') (27) 

aj 

\Jext, we consider the evaluation of the total angular distribution (eq. (12)). Thus, 
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da\ 

d£l ) XOT 


2 n | J ^ d 2 b' ex P ■ ( b - b 0] < OpOrj 1 1 - JJ [l - r a j (b - s Q - Sj 

l a.7 

x | 1 - n [i - r+.(b' -s Q - Sj )]J I OpO T > ■ 


( 28 ) 


which reduces to 


S T0T = I£I7 A dVexp [ * q ■ (b ~ b ' )] { 1 - n t 1 - ^(b)] - n i 1 - M i( b o] 

^ a J Otj 

+ n I 1 - M °J (*>) - Kj ( b ') + r aj (b, b')] 1 (: 

otj J 


where 


T oj(b. b ') = J dr a drj p P (r a )pT (rj) T aj {b - s a - Sj ) r+(b' - s a - Sj) (30) 

Combining the first two terms on the right-hand side of equation (13) gives 


da PT \ _ da p \ da T \ 

d / IN dQ.) IN dSl y iN dQ J 


(31) 


IN 


which, using equations (19) and (20), yields 


da 

dnj 


i = 


IN 


ik 

2n 


J d 2 b d 2 b' exp [iq • (b - b')] { II I 1 - M a,( b ) - M aj( h ') + W b - b')] 

v aj 

-110- [> - <(>>')] ] (: 

aj ) 


Equation (32) may be approximated by 

i2 


da\ 

dn) 


) = 
IN 


ik 

2tt| 


J d 2 b d 2 b' exp [iq • (b - b')] exp |i [*(b) - X + (b')] } {exp [T(b, b')] - 1} 


We can introduce momentum space representations into equations (20), (23), and (30) to find 

X(b) = / d 2 qexp(iq b) FW(— q) G (1) (q)/ NN (q) 

A 2 r>A 


(33) 
find 

(34) 


A LA p 

f2( b , b') = ^ 27 rfc N j^) 2 J d2 Q d2 d exp(iq • b) exp(iq' • b') 

x /NN(q) /nn(-q') F< 1 >(-q)F< 1 >(-q / )G( 1 >(q + q') 


and 


T(b, b') = d 2 q' exp(iq • b) exp(iq' • b') 

x /nn(q) /nnC-qO F (1) (q + q') G (1) (q + q') 


(35) 


( 36 ) 
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where and G (1) are the projectile and target one-body form factors, respectively. Equa- 
tions (34)— (36) may be more convenient for the evaluation of these phases than previous 
expressions. 

Correlations and Inclusive Scattering 

The effects of short-range dynamical correlations and Pauli blocking in the nuclear wave 
function will be most pronounced in the inelastic distribution at small and medium momentum 
transfers. In order to include these effects in the inelastic scattering distributions, we consider 
the ECC model. Assuming correct or equivalent kinematics, the ECC can be considered 
the matrix representation of the Glauber amplitude. In the ECC the matrix of scattering 
amplitudes for all possible projectile-target transitions is given by (ref. 17) 




d 2 b exp(zq - b) {exp[ix(b)] - 1} 


where barred quantities represent matrices and the elements of x are written as 

< mpa T |x(b)| n P vp >= — V f d 2 q exp(iq • bJF^pnpC-q) G>jW(q) /nn(q) (38) 

27 rfc NN J 

Assuming that the off-diagonal terms in x are small compared with the diagonal terms, we 
separate x into diagonal ( Xd ) an d off-diagonal ( Xo ) terms as 

X(b) = Xd( b) + Xo(b) ( 39 ) 

We further assume that the nuclear density in the excited states is approximately the same as 
the ground state, such that the elements of the diagonal matrix xd a-re all taken as the elastic 
element, 

x(b) = [ d2 q ■ F ’ (1) (-q) G (1) (q) /NN(q)exp(iq • b) (40) 

27 tA:nn J 

To treat off-diagonal scattering, we expand f in powers of xo> 


/( q) = f d 2 b exp(iq-b) exp [ixj^b)] ^ 

J 771=1 


[iXo(b)] T 


The inclusive distribution for the projectile then follows as 


J d 2 b d 2 b' exp [iq • (b - b')] exp ji - X + (b')] j 


x [T S (b,b') + T D (b,b , ) + ...} 


where the single inelastic scattering terms are 


Ts(b,b') = < OpOp |x(b)| Opv T >< v T Op x + (b') OpOp > 


and the double inelastic scattering terms are 

T d { b, b') = ^ E < °P®T lx(b)i V-pnp >< HT n P lx(b)| Opvp > 

0 np — 0 

X ^2 E < u tOp x + (b')| UpTi'p >< n'pn'p x + (b')| OpOp > 

n'p= 0 
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Each term in the inelastic scattering expansion of equation (42) can be reduced through the 
use of closure to terms involving matrix elements of one-, two-, . . . , etc., body operators over 
the ground state and thus includes the effects of two or more particle correlations. The single 
inelastic scattering terms may be reduced to 

E T s ( b, b') = <0 T \< Op\ x(b)|Op >< Op |x + ( b 0| °P > I °T > 

— < O pOp|x(b)|0 pOp >< OpO T |x + (b')|OpOr > (45) 


which becomes (ref. 16) 


E T si b ’ b ') = 
0 


^ p ^ T - s f d 2 q d 2 q' exp(iq • b) exp(iq' • b') 

(27rfc N N) J 


x / NN (q) F( 1 ) (q) ^ (1) W) 

X [^(q + q') + (At - 1) G (2) (q-q') - At G^fa) G^fa')] 


(46) 


where G (2) is the two-particle form factor of the target. A common approximation (ref. 18) 
is to neglect the renormalization of the one-body density to be consistent with the two-body 
density such that 

G^(q,q') = G ( 1 ) (q) G (1) (q') [l - C(q + q')] (47) 

where C(q -f q^) is a correlation factor. We then find 

r 5 (b,b') = fl(b,b') - j- x(b) X + (b') + K( b,b') (48) 

LS T /0 T 

where ft is defined in equation (35) and K represents the correlation term. The second and 
third terms on the right-hand side of equation (48) represent corrections to the model given 
by equation (25). We note that the second term in equation (48) persists even if the IPM is 
assumed for the nuclear wave function. 

The higher-order terms in equation (42) quickly become intractable as we go past the 
single-scattering term. Particle production multiplicities are generally small (< 3) for light- 
to-medium nuclei, suggesting that a perturbative approach would be useful, especially with the 
reasonableness of the distorted wave Born approximation (DWBA). A more fruitful approach 
is to look for a summation of the higher-order terms in a simplified, although approximate, 
manner. In order to make such a summation, we consider the double-scattering terms in 
equation (42). Prom equation (44) we have 



E T d( b.b') 


= E 


Uj'^0 


< OpO T Ix(b) x(b) I Opu T > 


- < OpOp |x(b)l Op >< Op Ix(b)| Opup >j 
< vpOp x + (b') X + (bO| OpOp > 

< ^rOp|x + (b')|Op >< Op |x + (b')| OpOp > 


(49) 
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which becomes 


£ 7b(b,b') = l -< 0 P 0 T \l x(b) x(b)| 0 P X 0 P \x + (b') x + (b') 

-x(b) X(b)| OpOj >< OpOr|x + (b') X + ( b ') 

- x(b) x(b)|0 P >< Op|x + (b0|Or >< O r |x + (b') 

+ X(b) x(b)|OpO r X 0 P 0 T \x + (b')\0 T X O r |x + (b') 

- x(b)|Op X Op|x(b)|Op >< Op|x + (b') X + (b') 

+ x(b)|O r >< Orlx(b)|OpO T x OpO r |x + (b') x + (b') 

+ x(b)|O r X O r |x(b)|0 P X 0p\x + (b')\0 T >< 0 T |x + (b') 

- x(b)|Op >< OT|x(b)|OpOp >< OpOp|x + (b , )|Op > 

<O r |x + (b')||OpO T > (50) 

To further reduce these terms, we write 

£ T D (b > b') = i£T Di (b,b / ) (51) 

vt^O i = 1 

and write the general term as 

T D .(b,b') = (- * ) J d 2 q\ d 2 q 2 d 2 q 3 d 2 q 4 

x exp [i(qi + q 2 ) • b] exp [i(q 3 + q 4 ) • b'] 

x /nn(qi) /nn(Q2) /nn ( qa ) /nn^) 

x Fi(qu q 2 , q 3 , q.j) <?i(qi, q 2 , q 3 , q4) (52) 

where F* and G x are functions of the projectile and target internal structure, respectively. 

The projectile function F* is identical in each term of equation (50) and is given by 

Fi{ qi, q 2 , q 3 ,q 4 ) = £ < Op|exp(iqi ■ s ai ) exp(iq 2 • s Q2 )|Op > 

x < Op I exp(iq 3 • s a3 ) exp(iq 4 • s ft4 )|Op > (53) 

which is written 

F t (qi . . .q 4 ) = [>4p F^(qi + q 2 ) + A P (A P - 1) F^(qi,q 2 )j 

x j^Ap F^)(q 3 +q 4 ) + Ap(Ap — 1) F^ 2 ^(q 3 ,q 4 )j (54) 
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We note that in the IPM we have 


*i(qi • • • q 4 ) — [Ap fW( qi + q 2 ) + A P (A P - 1) Qi) F^( q 2 )] 

x [A P F«( q 3 + q 4 ) + A P (Ap - 1) F^qs) F^(q 4 )] (55) 

If the coherent approximation is made for the elastically scattered projectile, the summation 
over all projectile states in equation (44) is not considered, and then we have 

Fj (qi . . .q 4 ) = A% FW( qi) F^( q 2 ) F<»( q 3 ) F^( q 4 ) (56) 

Next, we list the target structure terms: 

Gi(qi . . .q 4 ) = - E <°T I exp(iqi • Sj 4 ) exp(iq 2 • s, 2 ) exp(iq 3 • s J3 ) exp(zq 4 ■ Sj 4 )\0 T > 

3 1-34 

(57a) 


G 2 (qi • • • q 4 ) = 51 <0 T | exp(zqi • Sj, ) exp(iq 2 • s j2 )|O r > 
j 1-34 

x < O t | exp(zq 3 • s j3 ) exp(zq 4 • s j4 )|O r > (57b) 

G 3 (qi • ■ • q 4 ) = 51 <°t\ exp(iqi • s,-,) exp(zq 2 • s J2 ) exp(iq 3 - s j3 )|O r > 

3 1 ■■■34 

x < 0 T |exp(zq 4 • s J 4 )|O t > (57c) 

G 4 (qi . . . q 4 ) = E < O t | exp(z'qi • Sj 4 ) exp(zq 2 • s J2 )|Or > 
h- 3 i 

x < O t | exp(iq 3 • s J3 )|Or >< 0 T | exp(iq 4 ■ Sj 4 )\0 T > (57d) 

G 5 (qi--.q4)= E < O r |exp(iqi -SjjJIOr > 

3 1-34 

x < O r | exp(iq 2 • Sj 2 ) exp(zq 3 ■ s J3 ) exp(iq 4 • s J 4 )|O t > (57e) 

G6(qi--.q4)= E < O r |exp(iqi -Sj,)|O r >< 0 T |exp(iq 2 • s, 2 ) \O t > 

3 1-34 

x < O r |exp(zq 3 • s J3 ) exp(zq 4 • s j4 )| 0 T > (57f) 

G?(qi . . • q 4 ) = E < °t\ exp(iqi • s ix )\0 T >< 0 T \ exp(iq 2 • s h ) exp(iq 3 ■ s j3 ) \0 T > 

J 1-J4 

x < 0 T I exp(zq 4 • Sj 4 )|O r > (57g) 

Gs(qi • • • q 4 ) = E < O r |exp(iqi •s, 1 )|Ot’ >< O r |exp(iq 2 s j2 ) \0 T > 

31-3 4 

x < 0 T I exp(iq 3 • s J3 )| 0 T >< 0 T \ exp(iq 4 • s J4 )\O t > (57h) 
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which become 


C?l (qi . • • q4) = ~ 1 )(At - 2 )(At - 3) G^(qi, q 2 , Q3> Q4) 

+ (A t - 1 )(A t - 2) [G (3) (qi + q2,q3,q4) + Permutations] 

+ {At - l) [c (2) (qi 4 q2> q3 + q 4 ) + Permutations] 

+ {At ~ 1) [G^(qi, q2 + q3 + qt) + Permutations] 

+ G (1) (qi +q 2 + q3 + q4)l ( 58a ) 


G 2 (qi • 

■ • q 4 ) = 4 j G^ 1 - (qi -1- q 2 ) + {At - 1) G^(qi,q 2 )] 



x ^G-^(q3 + q4) + {At — 1) G^(q3,q 4 ) | 

(58b) 

G 3 (qi • 

. . q 4 ) = A\ G (1) (q4){(^T-l) {At- 2 ) G {3) (qi,q2,q3) 

4- {A t - 1) [G^(qi 4- q2, q3) + Permutations] 



x G (1) (qi + q 2 + q3)| 

(58c) 

G.i(qi 

. . q 4 ) = —A 3 G^(q 3 ) G^(q 4 ) jG^(qi + q 2 ) + {At - 1) G^(qi, q 2 ) 

'cn 

00 

Gs(qi 

..q 4 ) = G 3 (q 4 ,q2.q3,qi) 

(58e) 

Ge(qi 

• • q 4 ) = G 4 (q 3 , q 4 , qi, q 2 ) 

(58f) 

G7(qi 

• • q 4 ) = G 4 (q 3 ,q 2 .qi,q 4 ) 

(58g) 

G 8 (qi 

. .q 4 ) = 4 GW(qi) G^{ q 2 ) G^qs) G^(q 4 ) 

(58h) 

where and G^ are the target three- and four-particle form factors, respectively. 



The reduction of the single inelastic scattering terms as given by equation (46) contains 
the two-particle form factor, whereas the double inelastic terms in equations (57) and (58) 
display the two-, three-, and four-particle form factors. Thus, correlation effects may lead to 
two-particle knockout for a single inelastic scattering on a target nucleon, and to three- or 
four-particle knockout for double inelastic scattering. 

Upon identification of the double-scattering terms, and assuming the coherent approxima- 
tion for the projectile, the approximation of equation (25) is seen to contain only a single term 
in equations (58). A model that retains the dominant contributions of two-particle correlations 
in the double and higher terms, while assuming the coherent approximation for the projectile, 
is to assume 


da F 
d fi 


) = 

ik 

' IN 

27 r 


J db db' exp [iq • (b - b')] exp {« jx(b) - X + ( b ')] } 
x (exp [Qc( b > b 0] - 1) 


(59) 


with 

fl c (b,b') = Y, T s(b> b ') ( 6 °) 
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( 61 ) 


Upon comparison we find, within the coherent approximation, 

b') = - (T D2 + T Di + T d g + T d$ ) 

The approximation of equation (59) will thus treat the double and higher inelastic scattering 
terms in an approximate manner but should be accurate if the inelastic scattering series 
converges quickly. A similar analysis for the complete inelastic distribution as defined in 
equation (13) could now be made using the coupled-channels approach, but it will be addressed 
in future work. 

Model Calculations 

We now consider the evaluation of the inelastic distributions discussed above. Ignoring spin 
effects, we use an isospin-averaged, two-body amplitude given by 

/nn(q) = a ^ a "^^ NN - exp(-Bq 2 /2) (62) 

where the energy-dependent parameters a , B , and a are listed in table I. 

Table I. Parameters for Nucleon-Nucleon Scattering Amplitude 


Reaction 

cr, fm 2 

B, fm 2 

a 

a - a at 642A MeV .... 

3.93 

0.13 

-0.39 

a— 12 C at 3.64A GeV .... 

4.2 

.28 

-.43 

p- 16 0 at 1A GeV 

4.3 

.26 

-.22 

a — Ax at 1A GeV 

4.3 

.26 

-.23 


For the projectile, we use a one-body form factor 

F (1) (q) = exp(-R 2 p q 2 /4) (63) 

where Rp is the matter radius of the projectile. For the target one-body form factor, we use 
the harmonic- well form (ref. 19) 

G (1) (q) = (1 - C T q 2 ) exp(-i4<? 2 /4) (64) 

where R p is the matter radius of the target and 


c T = 7r *f 

4(1 + t) 


(65) 


with values of pp from reference 19. 

Correlation effects are included in the two-particle density through the approximation 
(ref. 18) 

p (2) (x,y) m p (1) (x) p (1) (y) [l - exp [-(x - y) 2 /2^]] (66) 

where £q, an effective correlation length, is equal to 0.7 fm. With the inputs of equations (62)- 
(66), the functions x(b)> ^c(b, b'), and Y(b, b') are evaluated in analytic form. For 
comparison with experimental results, the inclusive invariant distribution is written as 


d(J P \ ^ 7 r 

An ~* 2 " An 


li 



with 


( 68 ) 


t = -q 2 

In figure 1 we show comparisons with experimental results from reference 11 for a a 
scattering at 642 A MeV. The solid line denotes the correlation model of equation (59), the 
dashed line denotes the IPM of equation (24), and the dotted line denotes the IPM of 
equation (25). The correlation model produces good agreement over the region of momentum 
transfers studied. The IPM results overestimate the data, thus indicating the importance 
of correlation effects in this reaction. Inclusion of more detailed two-body densities in the 
calculation should improve the predictions. 

In figure 2 we show results for inclusive inelastic ot —^ 2 C scattering at 3.64A GeV. The 
correlation model shows a slight decrease in magnitude in comparison with the IPM results. 
Also apparent is a slight dip at small values of t. The increase in elastic coupling for increasing 
target or projectile mass number should mask correlation effects, a result that is seen to be 
true upon comparison with the data in figures 1 and 2. In figure 2 the approximation of 
equation (25) is seen to adequately represent the more exact IPM results of equation (24). 
Experimental results (ref. 20) for total inclusive a- 12 C scattering at 3.64A GeV are shown (see 
data points) in figure 3. The dashed line denotes the inelastic results of equation (25); the 
dotted line, the elastic contribution calculated using the coherent approximation described in 
reference 17; and the solid line, the sum of inelastic and elastic contributions. Agreement with 
the data is fair, where calculations underestimate the data at larger values of t. Correlation 
effects in elastic scattering have been shown to increase the cross section in this region by a 
substantial amount (refs. 16 and 17) so that use of a second-order elastic scattering model 
should lead to improved agreement. The dominance of elastic scattering at small values of t , 
as seen in figure 3, indicates that the model of equation (29) is sufficient when total scattering 
distributions are considered. 

In figures 4 and 5 we show results for inclusive inelastic p " O scattering at 1A GeV. 
In figure 4 the correlation model is seen to produce a substantial decrease in cross section in 
comparison with the IPM results. A much smaller decrease is seen for a 16 0 results in figure^. 

In figures 6 and 7 we show results for complete inelastic scattering in a a and a- C 
collisions at 1A GeV. The solid line denotes the exact IPM of equation (32); and the dashed 
line, the approximation of equation (33). This approximate form is noted to be inadequate at 
small momentum transfers. The complete inelastic distribution may provide a useful study of 
correlation effects at small values of q for identical projectile-target systems where the ground- 
state, two-particle form factor will play a dominant role in the leading-order terms. 

Finally, we note that the results of this paper require large computational times (2—3 hours 
per figure on a MICRO- VAX TS 05) because of the four-dimensional integrations over the two 
required impact-parameter points. Computational times will be substantially reduced when 
equation (29) is accurate. Here, an approximate reduction of the two-point impact-parameter 
integration may be used. 

Concluding Remarks 

The inclusive inelastic scattering of heavy ions is discussed in terms of the high-energy optical 
model using ground-state sum rules. The multiple-scattering structure of . these reactions is 
developed in terms of uncorrelated and correlated nuclear wave functions. Comparisons between 
the models are made and experimental data arc introduced using simple models for ground- 
state, one-, and two-body densities. Results indicate that correlation effects may be important 
only for proton projectiles and light-ion scattering. Improved physical inputs and numerical 
evaluation techniques are necessary for further applications. Approximate laboratory energy 
spectra for inelastically scattered projectiles can be derived from the invariant momentum 
transfer distributions discussed. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
August 24, 1990 
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Figure 2. Inclusive inelastic a- l2 C scattering at 3.64A GeV. 
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Figure 3. Total inclusive a- l2 C scattering at 3.64A GeV. 
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Figure 6. Complete inelastic a- a scattering at 1A GeV. 



Figure 7. Complete inelastic a 12 C scattering at 1A GeV. 
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